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NOIiENCLATORE 

■     '.  .      ■      -J 


>    .y, 


A 

2 
flow  area,  ft 

C 

sound  speed,  ft/sec 

^P 

constant  pressure  specific  heat,  ft  Ibf/slug  ®R 

D 

diameter  of  pipe,  ft 

F 

impulse  function,  Ibf 

f 

friction  coefficient 

g 

acceleration  of  gravity,  ft/sec^ 

k 

specific  heat  ratio 

L 

length  of  pipe,  ft 

M 

Mach  Number 

P 

pressure,  Ibf /ft  abs 

Q 

heat,  ft  Ibf/slug 

% 

parameter  defined  on  page  12 

R 

gas  constant,  ft  Ibf/slug  ^ 

T 

temperature,  ®R 

V 

specific  volume,  ft^/slug 

Z 

elevation,  ft 

P 

density,  slug/ft^ 

^w 

trail  shearing  stress,  Ibf /ft 

1* 

condition  at  Mach  Number  unity 

)>iit 

condition  at  Mach  Number  equal  to  1/  ^k 

'o 

stagnation  condition 

^h 

horizontal  property 

'v 

vertical  property 

f(  ) 

functions  defined  on  page  7 

!▼ 


■'  > 


g(  )   functions  defined  on  page  13 

h(  )   functions  defined  on  page  18  to  19 

j(  )   functions  defined  on  pages  23  to  2k 


.■*  f 


INTRODUCTIOH 

From  a  one-dimensional  point  of  view,  the  three  most  common  factors 
tending  to  produce  continuous  changes  in  the  state  of  a  flowing  stream  for 
horizontal  flow  are:  (1)  changes  in  cross-sectional  area  (2)  external  heat 
exchange,  and  (3)  wall  friction.  Furthermore,  if  the  gas  flows  along  a 
vertical  pipe,  the  effects  of  change  in  elevation  must  be  considered. 

Foxir  basic  cases  are  considered  in  this  report.  The  first  case  is 
isentropic  flow.  In  horizontal  isentropic  flow,  the  area  change  is  the  main 
factor  that  causes  the  change  in  pressure,  temperature  and  other  properties. 
In  vertical  isentropic  floxi  with  constant  cross-sectional  area,  the  main 
factor  that  causes  the  change  in  properties  is  the  change  in  elevation.  These 
two  flovra  have  been  correlated  by  a  functional  relation  between  area  change 
for  horizontal  flow  and  elevation  change  for  vertical  flow.  The  other  three 
cases  studied  are:  (1)  reversible  diabatic  flow  with  constant  cross-sectional 
area  for  both  the  horizontal  and  vertical  flows,  (2)  reversible  isothermal 
flow  Xvfith  variable  area  in  the  case  of  horizontal  flow  and  constant  area  in 
the  case  of  vertical  flow,  (3)  irreversible  isothermal  flow  in  ■lAich  the  area 
is  constant  in  both  the  horizontal  and  vertical  case. 

The  results  of  the  correlations  are  presented  on  graphs  in  irfiich  the 
flow  parameters  that  are  used  in  the  correlation  are  the  independent  variables 
and  other  quantities,  such  as  Mach  Number,  pressure,  etc.,  are  dependent 
variables. 

In  all  the  cases  studied  the  flow  is  subsonic »  and  in  the  vertical  cases 
the  flow  is  upward. 


CORREUTION  OF  VERTICAL  AND  HORIZONTAL  ISENTROPIC  FLOWS 

In  reversible,  isentropic,  horizontal  flow  it  is  the  change  in  cross- 
sectional  area  that  physically  causes  the  change  in  Kach  Number,  temperature, 
pressure  and  other  properties.  In  reversible,  isentropic,  vertical  flow, 
with  constant  cross-sectional  area,  it  is  the  change  in  elevation  that 
physically  causes  the  change  in  Mach  Number,  temperature,  etc,.  As  area 
change  is  the  fundamental  entity  that  causes  property  changes  in  horizontal 
flow,  and  elevation  change  is  the  fundamental  entity  that  causes  property 
changes  in  the  vertical  flow,  these  two  fundamental  entities  have  been  chosen 
as  the  basis  for  the  correlation  of  these  t\io  types  of  isentropic  flow.  As 
developed  later  this  correlation  is  obtained  through  the  equation 

m  k+1 

Correlation  of  other  properties  are  made  in  a  similar  manner.  For 

example,  the  correlation  between  (^)v  and  (~)  is  produced  from  the 

following  relations 

2(l-k) 
k+1 

1-k 

(1.)  =|l£ii  (II)     ^i^^ 

^T*V      Lk-1  ^T''h  '  k-lj 


and 


k+1 
4*''h  "^  k-rL4*V         -k+lj 


Isentropic  Horizontal  Flow 


The  relations  for  this  type  of  flow  are  t^ell  established  and  may  be 
found  in  any  standard  text  on  gas  dynamics.  The  eqxiations  used  in  this 
report  are  from  pages  83  to  87  of  reference  1. 


k+l 


k  +1 


4*^h 


2  + 


(k-l)l^ 


(1-1) 
(1-2) 


(p*)h 


(F^h 


=  L 


k  +  1   FI 


2  + 


(k-l)l>? 


J'- 


1  +  kl'i^ ^ 

ML2(k+l)(l+^M^)J 


(1-3) 


(1-^) 


Isentropic  Vertical  Flor^ 

Basic  Governing  Equations   For  isentropic  flow  in  a  vertical  pipe 
with  constant  cross-sectional  area,  the  basic  governing  equations  can  be 
written  in  the  following  forms: 

Energy  Equation 

CpdT  +  VdV  +  gdZ  «  0  (1-5) 

Momentim  Equation 
dP 


I 


Z  +^Z 

V  +  ^iV  I  M+^M 

1 


+  VdV  +  gdZ  =  0 


Continuity  Equation 


P   2  v2 


(1-6) 


(1-7) 


k 


Perfect  Gas  Relation 
P    p    T 


(1-8) 


(1-9) 


Definition  of  Mach  Nvunber 
d>f  ^  dvf   dT 

Working  Formulas   From  the  above  relations,  expressions  for  pressure, 
temperature  and  other  properties  in  tenos  of  Mach  Number  (M)  can  be  obtained. 
Combining  (1-5)  and  (1-6) 

C  dT  =  ~  =  C^  ^  =  kRT  ^ 
P     P      P       P 


v&ere  C  =  kR/k-1,  and 

f  =  ^  (k-1) 
T   p  ^ 

Substituting  (1-6) ,  (1-7)  into  (1-5)  and  eliminating  ^, 


,k+lx  dT 


(1-10) 


(1-11) 


Integrating  betueen  an  arbitrary  section  and  the  critical  section 


In 


M^ 


=  -OanT 


H 


k+1 

From  isentropic  relations 


T 


(1-12) 


(4)  =  (±)  k-i  _  ,i.vk+i 

From  the  definition  of  the  impulse  function 

/F_x  _  /_P  X  l-Hdl^ 
^j.*;^-  ^P*^v  1+k 


(1-13) 


(1-1^) 


Cosibining  (1-14)  and  (1-13) 

2k 


Similarly 


,  F  V    fUk+1   /1-HcM^N  (1-15) 


T  1  +  ^^  1^  ■ 

r-2^  a  fX.^  c         ^       )  (1-16) 


From  (1-12)  and  (l-l6) 

(1- 


T      2(l-k) 
(^)^  =  M      vi:j^  .  g^i . 


p         1  +  ^~1  M^  -iL. 

f_Ov       ,2.^    ( ?-.^)^"^  (1-18) 

o  1  +"2" 

PVom  (1-13)  and  (1-18) 

2k  k 

'  (S)  =  (J)^*"^  (4;  +  fer  ^i^)^"^  (1-19) 

^p*'v   M'    1-Hc   k+1   ' 
o 

j 
Relation  between  Z  and  M   Writing  the  energy  equation  between  any 

section  and  the  critical  section 

2   2.  ■ 

Cp(T*  -  T)  +  ^*  ^-  ^  +  gAZ*  =  0 

Dividing  through  by  CpT* 

/".       T  N    .    /k-lN  V*^  -  V^   .  k-1  AZ*  „  ^  n 


2 


(1-20) 


(1-21) 


From  (1-20).   (1-21)  and  (1-12) 

Correlation  of  Vertical  and  Plorizontal  Flows 

Rearranging  (1-22) 

2(l-k) 

1  +  i^  -S-  ^z*  =  (M)  ^-"^      (jIj)  (1  +  ^  m2)  (1-23) 

C* 

k+1  k+1 

The  right  side  of  (1-2*4-)  is  equal  to  the  area  ratio  (■^)  for  horizontal 
isentropic  flow,  thus 

k+l 
[l  +  ^  ^  tZ-f^  =  i^)  (1-25) 

This  relation  is  the  equation  that  is  used  to  correlate  the  vertical  and 
horizontal  isentropic  flovra.  Similarly,  from  (1-15) 

2k  1, 

and,  from  (1-4) 

,F_.   1  +  kl^ 

^^*^^"M[2(k4.1)(l+iS^M2)Ji 

Rearranging 

1  +  kM^  =  (^)j^  l(2(kfl)(l  +  ^  m2)J^  (1-26) 

From  (1-15)  and  (1-26) 

rX^  =  rl-^  L2  +  (k-l)I>I^/  ,    . 

(k+l)^(M)^^ 


From  (1-12)  and  (1-2) 

1-k 

/  T  V         I  k+l/_Ti>     -  _2_  P^  (1-28) 

^"t*  V  °  Lk-1^  T  ^h      k+lJ 

From  (1-13)  and  (1-3) 

Fron  (1-2).  (1-12)  and  (1-17) 

(-^)^  =  (-If)/-?  )b  (1-30) 

o 

From  (1-3),  (1-13)  and  (1-19) 

(^),  =  (-^),(-2),  a-31) 

From  equations  (1-1)  and  (1-22)  it  is  seen  that 

Bierefore  '  > 

(^)^  =  (-|j),f3L(-i^)hJ  =  (-li^vii  *  ^  sfe^^'J 

IAket.rf.se,  fron  equations  (1-2)  and  (1-22) 

From  equations  (1-3)  and  (1-22)  and  T-P  relation 

Similarly 


The  results  of  the  correlation  have  been  plotted  for  subsonic  upward  flow  and 
are  shoicn  in  Fig.  (2). 
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FJQ.    2 


5iO 


k-t-i 


-      /_A^)    _;.+  2(^^-1)    9      a2*1^<^^ 
I  A*/h  — l*^      ki-i      ^*  ^  J 

CORRELATION   OP   HORIZONTAL  $  VERTICAL 
ISEKITROPIC   FLOV/S  .__   . 


CORRELATION  OF  VERTICAL  AKD  HORIZONTAL  REVERSIBLE 
DIABATIC  FLOWS  IN  CONSTANT-AREA  PIPE 

Reversible  Diabatic  Horizontal  Flcn^  in  Constant-Area  Pipe 

The  relations  for  this  type  of  flow  are  well  established  and  may  be 
found  in  standard  texts  on  gas  dynamics.  The  equations  used  here  are  from 
pages  19^  to  196  of  reference  1. 

(X)  ^M^  (2-1) 

^T*^^   (Hrf)2 

T^    2(k+l)M2(l  +  ^  M^)  ^   ^ 

i^).-^^  C2.3) 


^P*'h 


1  +kl^ 

k 


(p|)h-^Li^-^M«^r  ^^-'^> 

(L\  =  -^  ^  (2-5) 

In  this  case,  it  is  the  heat  flow  that  causes  the  change  in  properties 
and  thus  it  is  the  heat  flow  that  is  taken  as  the  independent  variable. 

Reversible  Diabatic  Vertical  Flow  in  Constant-Area  Pipe  . 

Mien  comparing  reversible,  vertical  flow  in  constant-area  pipes  for 
the  two  cases  of  diabatic  flow  and  isentropic  flow  it  is  to  be  noted  that 
the  momentum  and  continuity  equations  are  the  same,  viiile  the  energy  equation 
differs  by  the  heat  flow  term.  In  the  isentropic  case  it  is  the  change 
in  elevation  that  causes  changes  in  properties.  In  the  vertical  diabatic 


10 


case  it  is  the  change  in  elevation  and  the  heat  flow  that  causes  the  change 

in  properties. 

In  reversible,  diabatic,  horizontal  flow  (Rayleigh  Line)  the  relation 
Pv^  equal  to  a  constant  is  not  the  polytropic  case  in  that  n  is  not  a 
constant.*  3Ji  order  to  get  a  solution  for  the  reversible,  diabatic  flow 
the  author  found  it  necessary  to  assume  that  n  is  a  constant  in  Pt"  equal 
to  a  constant.  Therefore  the  correlating  factor  in  vertical  flow  is  a 
function  of  heat  flow,  elevation  change  and  the  value  of  n. 
Basic  Governing  Equations 

Energy  Equation 

dQ  =  CpdT  +  VdV  +  gdZ        (2-6) 
Momenturi  Equation 


7  +  ^Z 
Y+AV  I  M+iM 

I 1 -: 


-f 


I 


Mk 


e,  V,T.  IW.Z. 

P16.  3. 

Definition  of  Mach  Number 


~  +  VdV  +  gdZ  =  0 
P 

Continuity  Equation 

Perfect  Gas  Relation 
dP   dp  ^  dT 

Tex  ■  *■  +  "ST* 
0  ^  T 


dK 


?"      dV^ 


dT 
T 


M^   V^ 
P  -  V  Relation 

Pv^  =  constant 


(2-7) 


(2-8) 


(2-9) 


(2-10) 


(2-11) 


*  It  can  be  shown  that  n  varies  frcan  zei^s  to  k  as  M  goes  from  zero  to  unity. 


■ye^-.  :.!"■ 
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^forking  Fortmilas   From  the  above  equations  by  assuming  n  a  constant, 
the  following  useful  formulas  can  be  obtained.* 

QJ«CpT*(l-X)-^RT*(l-|r)  (2-1^) 

,T.      5 (2-15) 

X^  =X4     .  (2-17) 


(t*^v  ~  (t*^v^  t  .  k-1  ^ 
o        ■>-  ■'■"T" 

Correlation  of  Vertical  and  Horizontal  FIotts 
Integrating  (2-6)  between  an  arbitrary  section  and  the  critical  section 
Q*  -  gAz*  »  C  (T*  -  T)  +  ^  [l  -  (^)^J  (2-22) 

Dividing  through  by  C*2 

From  (2-16) 

From  (2-15) 

*  Derivations  are  in  APPendix  I.  .        • 
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^T^K  =  ^  -  — : — n 


CpT*  -  ^  RT* 


As  C    =  ^  and  C*^  =  kRT*.   (2-15)  becomes 

{±\    -  1       {'^\    I      k(k-l)(n-l)      I 
(^)^  -  1  -  i.-^^  Lk(n-l)  -  n(k-l)J 


C* 
Let 


-ys^^^j 


(^)^  =  1  -  (^)^  N  (2-25) 

Combining  (2-24)  and  (2-25) 

Combining  (2-25),  (2-26)  and  (2-23) 

c*         c*  c* 


=  (^  -  Fx^'v  -  -^ 


C*^  V      2 


I  Q^  ■  g^2*  ,  >^  -  l-j.C*2   ^2-  (k-l)M2  . 

L   ^^2    +    2^J  WV  °   2(k-l)  ^2"^^^ 


C* 


From  the  horizontal  relation 
Combining  (2-1),   (2-5)  and  (2-28) 

Simplifying 

(^)     =     ("^DO^-f  (2.30) 
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Comparing  (2-2?)  and  (2-30) 


CX)    -I^Llf^+l^k^)  {^^^){^-^f  (2-31) 


From  (2-30)  and  (2-31) 

M  =  gi(Qj)  =  ggCQJ.AZ*) 
From  (2-15)  and  (2-31) 

(■§p\  =  g3(Qg)  =  g4(Qj.AZ*) 

From  (2-15),  (2-17)  and  (2-31) 

(^)^  =  g5(Qg)  =  g6(Q*.AZ*) 


From  (2-15),  (2-18)  and  (2-31) 

(||)^  =  g/Qj)  =  ggCOf  .^Z*) 
From  (2-17),  (2-19)  and  (2-31) 

(^)^  =  g9(QJ)  =  gioCQj.^Z*) 
Since  (|y)^,  (|5)j^,  (^)j^,  (^)j^  and  M  are  functions  of  (-^1)^,  and 

(|t)^.  (pp)^,  (~)^,  (^)^  and  M  are  functions  of  (-g)^  and  (-2^  AZ*).  the 

two  flows  are  correlated  through  equation  (2-31).  !I5ie  results  of  this 
correlation  are  shown  on  Fig.  (4). 
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q7s\ 

-A.  M 

\'^ 

^^-"""(P/R^-^ 

2.0 

0.50 

\ 

/ 

X 

^^ 

"^n^Ti^;^-^ 

- 

0.25 

^""^v^,.,^^ 

^'^ 

\      ■ 

I.O 

""\mT«),_ 
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\ 

^ 
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■c«V     1 


C*2         "^  2 


2.0 


3.0 


^C*^;^       I         c*2        '^         2     J    \Q*N/V      (n.kM*)*i:2-(k-l)M»] 
Fl<5.  4-    CORKELATIOKi  Ol^    HOKIZIOKITAL     ^    V&RTICAL 
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CORRELATION  OF  HORIZONTAL  AND  VERTICAL 
REVERSIBLE  ISOTHER>IAL  FLOWS 


In  this  part,  reversible  isothermal  cases  in  both  horizontal  and 
vertical  flows  are  considered. 

Reversible  Isothermal  Flow  in  Horizontal  Pipe  Td.th 
Changing  Cross- Sectional  Area 

The  properties  of  gas  flow  in  horizontal  pipe  with  changing  cross- 
sectional  area  are  effected  by  the  area  ratio.  For  the  isothermal  case, 
heat  transfer  occurs,  but  thie  flow  of  heat,  in  order  to  maintain  constant 
temperature,  is  dependent  on  this  area  change.  In  this  case,  the  governing 
equations  become 


Basic  Governing  Equations 


Pva.  51 


Energy  Equation 

dQ  =  VdV 
Momentum  Equation 


^  +  VdV  =  0 


Continuity  Equation 

da+dA^dV^O 
p  ^  A        V 

Perfect  Gas  Relation 

dP  _  dp 

Definition  of  Mach  Number 

dM^  _  dV^ 
m2"  v2 


(3-1) 


(3-2) 


(3-3) 


(3-^) 


(3-5) 
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Iforklng  Formulas  From  the  above  equations,  useful  working  fonmilas  for 
horizontal  isothermal  flow  can  be  obtained.  ♦* 

(JL)^  =  JkM       ■  (3-11) 


v*t  ^ 


(3^) 


1  -  kM^  k 

o 

o 


<5*  =  f^  (1  -  kM^)  (3-17) 

(-^),  -  i^  (3-19) 

F**  ^   2  JkiM 


Reversible  Isothermal  Flow  in  Vertical  Pipe  ^Ath 
Constant  Cross-Sectional  Area 

In  this  case  the  main  factor  causing  change  in  properties  is  the 
change  of  elevation.  In  order  to  maintain  constant  temperature  the  flow 
of  heat  must  be  adjusted  to  the  change  in  elevation.  Therefore  it  is  the 
elevation  change  that  is  the  correlating  factor  in  this  case. 

Basic  Governing  Equations   For  a  reversible  isothermal  vertical  flour 
with  constant  cross-sectional  area,  the  governing  equations  become 

**  Derivations  are  in  Appendix  H. 
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y\\?     z  +  ^z 

V  +  AV  I  M+^M 


L 


,aQ 


Energy  Bjuation 

dQ  =  VdV  +  gdZ 
Kbmentum  B^iiation 

~  +  VdV  +  gdZ  =  0 
P 

Contimiity  Bjuation 

P   V 

Perfect  Gas  Relation 

P        P 


(>20) 


(3-21) 


(3-22) 


(3-23) 


Definition  of  Mach  Nuniber 

dll^      dV^ 
1?"    V^ 


(>2i^) 


Vforking  Formulas   From  the  above  equations,  useful  working  formulas 
are  obtained;* 


(-7?)^  =  fk  M 


V* 


t^v 


(3-30) 


5Az**4ln(^)-i(l-kM2) 


C*' 


/_P_N        1 

P**'^  ^   /k  M 


(^). 


bL^a-¥^^)7 


k 
:-l 


p*t'v   ^  jj  L3k. 

(r:^)v  =  31^  CI  +  -r  r) 


o 


(3-31) 


(>32) 


(3-33) 


(3-3^) 


*  Derivations  are  in  Appendix  U. 
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Correlation  of  Vertical  and  Horizontal  FIot-ts 
From  (3-13) 

t   2   ^ 


Substituting  into  (3-31) 

(j<,t  A*  A* 

From  (>11)  and  (3-13) 

V  V     (   2   '  A*^ 


From  (3-13)  and  (3-12) 


'^)  =  (-^ 

A* 

From  (3-13)  and  (y-V^) 


(— ^h  =  <-^>h  '^  ^  (3-38) 

A 


(|),=  (^),-^[^a.i^^)r  039) 

From  (3-13) 

^  =  ^1  L^-t^hJ  (3-^0) 

A* 

From  (3-29) 


-<^f-    "'!'-    ' 


2$ 


From  (3-31)  and  (3-32) 

C* 
!nierefore,  from  (>36)  and  (3-^) 

C* 
From  (3-36)  and  (3-^2) 

From  (>36),  (3-33)  and  (3-32) 

o        C* 

P       T 

Since  (-~)^,  (-~)^.  (-~)^.  and  M  are  functions  of  (-^).  .  and 

p*t  n   p:st  a   m*t  n  j^#^  n 

P      T 
(-L.)  ,  (-£-)  ,  (-2-)  ,  and  M  are  functions  of  (-%  2**),  the  two  flows 

p#t  V   pi«t  V   j«t  V  Q*C 

are  correlated  through  equation  (3-3^).  The  results  of  this  correlation 
are  shown  CEuFig.  (?)• 
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CORREUTION  OF  VERTICAL  AKD  HORIZONTAL  IRREVERSIBLE       /  ' 

ISOTHERMAL  FLOWS 

For  irreversible  case,  wall  friction  is  the  main  factor  Tdiich  causes 
properties  change  in  both  horizontal  and  vertical  flows,  the  heat  flow  is 
dependent  on  the  wall  friction.  Therefore,  the  correlation  of  these  two 
floira  is  based  on  this  fact. 

Irreversible  Isothermal  Flor-r  in  Horizontal  Pipe 

For  this  case,  relations  are  well  established  and  can  be  found  in 
any  standard  text  of  gas  dynamics.  The  relations  used  here  are  from  pages 
180  to  181  of  reference  1. 

o 
<;g*  =  ^  (1  -  ktj^)  (4-6) 

Irreversible  Isothermal  Flox^-  in  Vertical  Pipe 
In  case  of  combining  effects  of  elevation  change,  heat  transfer  and 
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friction,  the  governing  equations  become* 
Basic  Governing:  Equations 


X 


y\^?    z+az 

V4-AV  i  M-I-4M 


r 


L 


1 


Momentum  Equation 


ciQ 


Energy  Equation 

dQ  =  VdV  +  gdZ 
Continuity  Bijuation 

p        V       " 
Perfect  Gas  Relation 

P        P 

Definition  of  Mach  Number 

dM^  _  dV^ 
M^   T" 


(4-7) 


(^W8) 


(^-9) 


(^10) 


MP  +  pAgdZ  +  pAVdV  +  T^nDdZ  =  0 
Vforkins;  Formulas   From  the  above  equations,  the  expressions 
various  properties  can  be  obtained. 

^f.2. 


(4-11) 


for 


AZ*  = ^T>  "■  In rs -  -TiT-  In  Idl 


2g+fEI  '  8» 


4f  .  _^t   2RTf  -   /V*^x''  ^  ,2RTf  ^  , .  , 
D"  ^2*  =  -^  In  (-v-)^  +  (r^  +  1)  In 


2g  ^'i^ 
2g  +  ^V*t^ 


4jf   p*t  2 
^  2g  +  ^RT 


gD 


P  V 


Q^*  =  ^  (1  -  kl^)  -  gAZ** 


(4-19) 


(4-21) 


(4-23) 


(4-24) 


*  Derivations  are  in  Appendix  HI. 


>-ij?f;^:  ,?'■'•■ 
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Correlation  of  Vertical  and  Horizontal  Flo^ra 

In  this  case,  —  L  is  the  main  factor  in  horizontal  flow  -which  causes 
the  properties  change.  In  vertical  flow,  it  is  -^  Z  ^ahldh  causes  the 
properties  change.  Through  those  two  factors,  the  correlation  of  the 
vertical  and  horizontal  flows  are  established. 

From  (4-19) 


In  klf  =  ( ^)  In     2^   -  ^  T  ^2*        (4-25) 


2g  +^T 

From  (4-5) 

4f  ,     1  -  k>^ 


From  (4-2)  and  (4-5) 


(^'h  -  3i  (t  « 


From  (4-5) 

From  (4-3)  and  (4-5) 

o 

From  (4-4)  and  (4-5) 

,^o  V    .  ,4f  _. 
o 


(4-26) 


Combining  (^25)  and  (4-26) 

TW"-;i2--*'(^-^2|g>^  o^^<.  -2RgT^Z*    (4-27) 


2k 


From  (^6)  and  (^5) 


Qg*  =  is  Ct  I) 


From  (^19) 


M-;i6(T^> 


From  (^21) 


From  (4-23) 

From  (4-2it) 

,4f 


Q?*  =  ^5  <T  ^) 


Since  (-—)..  M.  (-—)..  (-%),,,  Q  are  functions  of  i^  L,  and  (—:)„, 

^n  >.  0  0 

P  T 

M,  (-tt)  ,  (-~r)  »  Q  are  functions  of  ~  Z,  this  correlation  between  these 

two  flows  is  established  through  equation  (4-2?).  The  results  of  this 
correlation  are  shoim  on  Fig.  (9). 
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APPENDICES 
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APPENDIX  I 
DERIVATION  OF  WDRKIKG  FORl-IULAS  FOR  REVERSIBLE,  DIABATIC  VERTICAL  FLOW 
In  this  case,  by  assuming  n**  in  ;Pv'^  =  constant  is  a  constant,  we  are 


From  (2-6)  and  (2-7) 


able  to  integrate  the  term  S£  in  the  momentum  equation. 


dQ  =  CpdT  -  ~ 


^  (2-12) 


1 


Pv»  =  const       v=  i^f  (2-13) 

Combining  (2-12)  and  (2-13) 

1 

dQ  =  CpdT  -  i^f  dP 

Integrating  between  an  arbitrary  section  and  the  critical  section 

1     n-1   n-1 
■  QJ  =  Cp(r*  -  T)  -  (?  ^(P*  ^  -  P  ^  ) 

■   =  Cp(T*  -  T)  -  ^  R(T*  -  T) 

«  CpT*(l  -  ^)  -  5^  RT*(1  -  ^)  (2-1^) 

From  (2-1^) 

C„T*  -  ~T  RT* 
P    n-x 


From  definition  of  I-Iach  Number 

1^ 
m2 


^T  N       /V  >2  i 

^Wv  =  ^F^v  «2 


**  Bxe  value  of  n  must  be  of  such  magnitude  that  there  is  heat  flow  into, 
the  flowing  gas.  This  assmaption  is  made  so  that  the  direction  of  heat 
flow  is  the  same  in  both  the  horizontal  and  vertical  cases. 


•^S--' 


^ 


p    n-1 


^  +  VdV  +  gdZ  =  0 
P 


Integrating 

—•  R(T*  -  T)  +  -""  I  ""  =  gAZ* 


V*^  -  V^ 


^  RT.(l  .  |r)  +  2^  L^  -  (^)^J  =  gW 


Dividing  through  by  C*^ 
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(^)  =m[i. \ J  (2-16)  • 

p    n-JL 

(?y)v  =  (p«)v(F^  =  ^F^J^ 

1  h 2L J^  (2-17) 

p    n-1 

(.f;)y  =  ( jy)y  (    k-1  ^ 
o        1  +-2- 

p    n-1 
(^)  =  (p*)v  (    k-1  ^ 


From  (2-7) 


j^(l-|,)4[l-(^)'j  =  ^M»  (2-20) 
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APPENDIX  n 
■  DERIVATION  OF  VWRKING  FORMULAS  FOR  HORIZONTAL  REVERSIBLE  ISOTHERMAL  FLOW 

In  this  case,  dQ  =  TdS  is  applicable,  and 

dS=.C  ^-R^ 
^   ^p  T     P 

dT  =  0 
dS--Rf 

-    ■  dQ  =  -RT^  (3-6) 

From  the  governing  equations  and  (3-6),  expressions  for  pressure,  temperature 
velocity,  etc,  can  be  derived  in  terms  of  M, 
Combining  (3-1)  and  (3-6) 

,  -  RT  ^  =  VdV  (3-7) 

Integrating  betxreen  an  arbitrary  section  and  the  critical  section** 

om  n       f^*'^\  V*^     -  V^ 

-  RT  In  i-y-)^  = ^ 

From  the  definition  of  Mach  Number 

V**  =  l>I*^A/kRT  (3-9) 

l-I**  =  1/  //k 

V**  =  aTrt 

•♦  Ihe  critical  section  for  reversible  isothermal  flow  occurs  at  the 
minimum  cross-sectional  area.  At  this  section  M  =  l/  //k,  (see  page  180 
and  problem  6-10  on  page  188  of  reference  l). 
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Coobining  (3-8)  and  (3-11)  ,,:^ 

/1-kM^v 

From  the  continuity  equation 
'  +   -fc      ^  2  ' 

^„£!:v!!^jg 1-  (3^3) 

From  the  definition  of  stagnation  pressure 

p  l+^M^  ^ 

.  ^;j^h  =  ^'^^h   ^1  ^.  k=l  1  ^ 
o  2  k 

.1-kM^N  k 

Froia  the  definition  of  stagnation  temperature 

o  2    k 

From  (>1) 

dQ  =  I  dV^ 


Integrating 


«£*  =  ^Li-<^''J  <>^« 
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From  (3-16)  and  (3-11) 

Q**  =  f  (1  -  kM^) 
From  the  definition  of  the  impulse  function 

(p^t^h  -  ^p«t^h^^.,t^h   2 
From  (3-12),  (3-13)  and  (3-18) 

.  F  V  ^  1  +kl-l^ 
F**  ^    kl-l(2) 


(3-17) 


(3-18) 


(3-19) 


DERIVATION  OF  VWRKING  FORI-IULAS  FOR  REVERSIBLE  ISOTHERMAL  VERTICAL  FLOW 


Fr<»i  (3-20)  and  (3-21) 

HO  -   ^  -   dP  P  _  _  j^«  dP 
dQ=-y=--jr---RT  p 

From  (3-25)  and  (3-20) 

-  RT  ^  =  VdV  +  gdZ 

From  (3-22)  and  (3-23) 

dP   dp  _   dV 
-y^   p  -  -  V 

ET  ^  =  VdV  +  gdZ 


(>25) 


(3-26) 


dV 


RT  ^  -  VdV  =  gdZ 


RT  f  -  I  dV^  =  gdZ 


(3-27) 


Integrating  between  an  arbitrary  section  and  the  critical  section 


RT  In  V 


V 
t 


-|(v2) 


Y*t 


=  g(Z**  -  Z) 


RT  In  (^)y  -  I  (V*"^  -  V  )  =  gAZ*^ 
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^  (^V  - 1  v*'a  -  -^)  =  g«'*  <>^> 


RT 

V* 


Dividing  through  by  C*^ 


^=  /kM  (3-30) 

Froa  (3-29)  and  (3-30) 

I  In  (--^)  -  ^  (1  -  kl'I^)  =  -^  AZ**  (3-31) 

From  the  continuity  and  the  perfect  gas  relation 

P*^  ^   //k  M 
From  the  definition  of  stagnation  pressure 

\^tW  -  ^p,tV  ^1  +  ki  1  ^ 

^o      ^    -^  ^  2  k 

k 

Similarly,  the  stagnation  temperature  is 

(""T)^  ~  (  t'v  (    k-1  1  ) 

Tjtv   ^tv  i+k_ll 

=  ^  (1  +  —^  M^)  (3-3^) 


?t'.;i 
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APPENDIX  III 


DERIVATION  OF  vroRKING  FORMULAS  FOR  VERTICAL  IRREVERSIBLE  ISOTHERMAL  FLOW 


2 

becomes 


Dividing  through  by  pA,  and  introducing  t^  =  £j-  f .  equation  (4-11) 


^  +  gdZ  +  VdV  +  ^(^)dZ  =  0  (^12) 

f  |  +  gdZ  +  VdV+^(^)dZ=0  (^13) 

From  the  perfect  gas  relation 

p 

£=  RT 
P 

(4-12)  becomes 

RT  ^  +  gdZ  + 1  dV^  +  ^  ^  dZ  =  0 

■  dz(g+^#)+RTf +  |dV^=0  (4.14) 

Rearranging 

-(RT^+|dV^) 
dZ  =      \.  \ (^15) 

From  (4-7),  (4-8)  and  (4-9) 

dP   dp  _  1  dV^  _   1  dl-I^  twjxf\ 

Combining  (4-15)  and  (4-l6) 

RT  dI-1^^  C^  ..2 

2  ^^      2  ^ 
dZ  = ?  o  ' 

,  .  4f  c'^r^ 


3«^ 


dZ  = 


RZ 


dI4 


2g^i^c¥    1.3^       2g+^c¥ 


dM^ 


(^17) 


Integrating  between  an  arbitrary  section  and  the  critical  section 


".**- 


AZ*t  =  ^ln 


Z    =- 


RT 


RT 
2g 


^l 


2g  + 


i^f 


M^ 


^-..^/'^ 


M 


2g+^cV 


D 


^        kM^  (2g  +  ^  RT) 


^1^ 


In 


2s^igcV 

2g+^RT 


4f    .„^t       RT  ^^  n 


2g^^c¥  2g-.iiic¥ 

kM^  (2g  +  ^  RT) 


+  ln 


-     .  4f  .2.  J2 

2RTf  4-  gP  ,„  '^^       D    ^  ^^         2RTf 


D 

5f 

2g  +-5-RT 

In  kl^ 


l/fk 


M 


(J^18) 


(^19) 


From  (4-17)  and  (4-10) 

RT 


dV^ 


^■iTTf^l^-^Tfv^ 


IV^ 


(4-20) 


Ihtegrating 


,,.t  =  |^(!il| 


^v^ 


-)C 


V*' 


2g  +  -^  V*^ 


) 


.2  2g+^V^ 


(4-21) 
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Fron  (4-16) 

Combining  Ot-ZZ)  and  (V19).  eltminating  M, 

«  ,„t .  aR«^  ^  ^"^;;^^'  -  ^  in  (^)^     (W3) 

^         2g  +  —■  RT       ^ 


From  (^7) 


dQ  =  I  dV^  +  gdZ 


Integrating 

Q«t  =  1|2  (1  -  kM^)  +  gAZ**  (^2i*) 


-?;•. 
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The  main  effort  of  this  report  is  to  find  out  the  relations  betwen 
gas  flows  in  horizontal  and  vertical  pipes. 

The  horizontal  relations  \iiXch  are  >rell  established  and  can  be  found 
in  any  standard  gas  dynamics  text.  The  vertical  relations  are  developed 

in  this  report. 

Four  cases  are  considered,  namely:  (1)  isentropic  flow  (2)  reversible 
diabatic  flow  (3)  reversible  isothermal  flo;^^  and  W   irreversible  isothennal 

flOVT. 

In  the  first  case,  change  in  cross-sectional  area  for  horizontal  flow 
is  the  main  factor  that  causes  the  change  in  pressure,  temperature  and 
other  properties.  In  isentropic,  constant-area  vertical  flow,  change  in 
elevation  is  the  main  factor  that  causes  the  change  in  properties.  Through 
these  Ujo  factors,  relations  between  these  two  flows  are  established. 

Likexase,  the  other  three  cases  are  correlated  in  a  simlar  manner.   . 


